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ABSTRACT 

The origin of dust in a galaxy is poorly understood. Recently, the surveys of the Large 
Magellanic Cloud (LMC) provide astrophysical laboratories for the dust studies. By 
a method of population synthesis, we investigate the contributions of dust produced 
by asymptotic giant branch (AGB) stars, common envelope (CE) ejecta and type II 
supernovae (SNe II) to the total dust budget in the LMC. Based on our models, the 
dust production rates (DPRs) of AGB stars in the LMC are between about 2.5 x 10 -5 
and 4.0 x lO _ 6 M 0 yr . The uncertainty mainly results from different models for 
the dust yields of AGB stars. The DPRs of CE ejecta are about 6.3 x 10~ 6 (The 
initial binary fraction is 50%). These results are within the large scatter of several 
observational estimates. AGB stars mainly produce carbon grains, which is consistent 
with the observations. Most of dust grains manufactured by CE ejecta are silicate and 
iron grains. The contributions of SNe II are very uncertain. Compared with SNe II 
without reverse shock, the DPRs of AGB stars and CE ejecta are negligible. However, 
if only 2 % of dust grains produced by SNe II can survive after reverse shock, the 
contributions of SNe II are very small. The total dust masses produced by AGB stars 
in the LMC are between 2.8 x 10 4 and 3.2 x 10 5 Mq, and those produced by CE ejecta 
are about 6.3 x 10 4 . They are much lower than the values estimated by observations. 
Therefore, there should be other dust sources in the LMC. 
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1 INTRODUCTION 

The interstellar medium (ISM) of a galaxy, one of the 
drivers of itself evolution, determines many characteristics 
of the stars forming in it. Dust is one of the important 
ingredients of the ISM, and plays a central role in the 
astrophysics of the ISM. However, the dust origin in the 
ISM is less understood. According to a popular point of 
view, the stellar winds of asymptotic giant branch (AGB) 
stars and type II supernovae (SNe II) ejecta hav e long b een 
considered the main site s of dust formatio n fe.g.. lGail et al.l 
120091 : [Dunne et ai1l2003l : iGall et al.l 1201 il l . In addition, red 
supergiants, luminous blue variables, Wolf-Ra yet stars and 


type l a supernovae can also p r oduce dust (e.g., MassevetjdJ 


20051: iMcDonald et alJ I2009I . l201ll : iBorkowski et al.l 20061 : 


Nozawa et al. 20 llll . However, compared with AGB stars 


and SNe II, th eir contributions to the dust production are 
negligible le.g.. IZhukovska et al JLOQsl ; IGall et al.ll201ll l 

Recently, as an ideal laboratory for dust formation and 
evolution, the Large Magellanic Cloud (LMC) was observed 
by a number of surveys whic h detected th e IR d u st emissions 
from stellar sources (e.g., ICioni et~al1 l2000l ; lEgan et~al] 
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200l]^ Meixner_et ah 20061 : iBlum et al Jl2()od : Ilia et al.ll2008l : 
Shimonishi et al.|[2013i l . Based on the observati ons from the 


Spitzer Space Telescope, iMatsuura et alJ (120091 1 carried out 
the measurement of the dust production rate (DPR) by 
stars for the entire LMC. They estimated that the DPR 
of carbon-rich (C-rich) AGB stars is ~ 4.3 x 10 _5 Mq yr _1 , 
much higher than that of oxygen-rich (O-rich) AGB stars 
(up to ~ 1.7 x 10 _5 Mq yr _1 ). Based on the best fit Grid 
of Red Supe r giant and Asymptotic Giant Branch Models, 
iRiebel et ali (|2012l l estimated that DPR by C-rich AGB 
stars is ~ 1.3—1.4 x 10~ 5 Mq yr _1 , which is two and a half 
times as much dust into the ISM as do O-rich AGB stars. 
Very recently, based on the evolutiona r y models o f stars 
through the AGB phase, iDell’Agli et al.l (|2014l . I2015T) inter¬ 
preted the Spitzer observations of extreme stars and AGB 
stars in the LMC. Theoretically, acc ording to the model for 
du st-m ass returned from AGB stars, IZhukovska fc Hennina 
J 2 OI 3 I) calculated the DPR of AGB stars in the LMC. They 
found that DPR by C-rich AGB stars is ~ 5.7 x 10~ 5 Mq 
yr _1 , and DPR by O-rich AGB stars is only ~ 1.3 x 10 -6 Mq 
yr _1 . Similarly, using the grids of dust yiel ds for different 
stellar masses and met a llicitie s calculated by Ventura j?^aL| 
(l2012fl. I Ventura et al.| 1 201^L_[DjAlrjsdenG_et ahl 1 2013lb 
and Ventura et al.l ( 2014 ), Schneider et al. ( 20141 1 found that 
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DPR by C-rich AGB stars is ~ 4.0 x 10 _5 A4'q yr _1 (~ 
3.0 x 10 _6 Mq yr _1 ) and DPR by O-rich AGB stars is 
~ 6.0 x 10 _6 Mo y r 1 (~ 6.0 X 10~ 6 Mq yr _1 ), respectively. 

The chemical composition of the dust mixture from 
AGB stars is determined by the carbon-to-oxygen element 
abundance ratio (C/O) in the stellar wind. In O-rich 
environment ( C/O < 1), dust produced is mainly silicate, 
while amorphous carbon dust is dominantly formed in 
C-rich environment ( C/O > 1). Based on the above 

observational and theoretical results, dust in the LMC 
should be mainly composed of amorphous carbon grains, 
and the silicate grains accou nt for a small part. However, 
IWeingartner & Drain j feOOll l estimated the total volume 
per H atom in the carbonaceous and silicate grain popu¬ 
lations in the Milk Way, LMC and SMC. The volumes per 
H atom of the carbonaceous and silicate grain in the LMC 
are ~ 0.4 x 10 -27 cm 3 H _1 and ~ 1.3 x 10 -27 cm 3 H _1 , 
respectively. Considering the ideal graphite density of 2.24 
g cm -3 and the silicate density of 3.5 g cm -3 , we estimate 
that the mass ratio of carbon to silicate grains in the LMC 
is about 1:4. 

Therefore, most of silicate grains may not be produced 
by AGB stars. Are they produced by SNe II ejecta? Un¬ 
fortunately, it is presently not defi nitely known how much 
dust is produced by SNe II ejecta. IZhukovska et al] (l2008l l 
simulated the dust returned by different sources in the solar 
neighbourhood. They found that the dust grains produced 
by SN ejecta mainly are carbonaceous grains, while the 
number of silicate grains is negligible. It is possible that 
most of silicate g r ains d o not originate from SN ejecta. 
iDraine &; Salpeten J 1 971)1) put forwar d that dust can grow 
in the ISM (also see IPrain e l2009l l . Unfortunately, the 
destruction and growth processes of dust in the ISM are 
poorly understood, and it is very difficult to observe them 
(e.g., IPwek el alj 2007). _ 

Recently. iLii et aD (l2013l ) suggested that the common 
envelope(CE) ejecta in close binary systems can provide 
a favourable environment for dust formation and growth. 
Very recently. I Nicholls et al.l ll2013l l showed that V1309 Sco 
had become dominated by mid-IR emission since eruption, 
which indicated the presence of a sign ificant amount o f 
dust in the circumstellar environment. I Zhu et ahl d2013l l 
considered that these dust grains around V1309 Sco were 
efficiently produced in its progenitor-binary-merger ejecta. 
The discovery of these dust grains offers an indirect evidence 
for the dust formation and growth in the CE ejecta. 
Therefore, CE ejecta is an efficient dust source. However, 
to our knowledge, there is no work considering CE ejecta as 
dust source on estimating dust budget in the LMC. 

In this work, we investigate the contributions of dust 
formed by CE ejecta to the interstellar dust, and give a new 
dust budget in the LMC. In §2 the model of the ejecta during 
binary merger is described. The main results and discussions 
are presented in §3. §4 gives conclusions. 


2 MODEL 

Before discussing dust formation and growth in CE ejecta, 
we give some descriptions for CE and CE ejecta evolution, 
and the chemical abundances of CE ejecta. 


2.1 Common Envelope Evolution 

As a result of dynamical timescale mass exchange in close 
binaries, CE evoluti on plays an ess e ntial role in binar y 
evolutions (see, e.g.. iPaczvnskil Il976l : llben fe Liviol 1 1 9931 ) . 
Unfortunately, based on observation data and theoretical 
models, CE evolution is poorly kno wn. In this work, we 
adopt the scheme by iHurlev et al.l (120021 ) to simulate CE 
evolut ion. All details can be seen in Section 2 of lHurlev et al.l 
(120021) . In most cases, CE evolution involves a giant star 
(donor) transferring matter to a normal or a degenerate star 
(gainer) on a dynamical timescale. If the orbital energy of 
the binary system is large enough, the whole envelope of the 
gainer is ejected as CE ejecta. 


2.2 The Evolutions of Gas Density and 
Temperature in CE Ejecta 

When CE ejecta begins to expand, its temperature and mass 
density star t to decrease. Acco rding to the classical theory of 
nucleation (iFeder et al.lll966l ). in most cases the formation 
of dust grains can occur in adiabatically expanding gas, 
and is mainly determined by temperature, mass density and 
chemistry of the gas phase. In order to simulate the dust 
formation and growth in CE ejecta, we must determine 
the evolutions of gas density and temperature. To our 
knowledge, there is no any comprehensive theoretical model 
or observational datum to describe the mass density and the 
temperature of CE ejecta. 

Considering the whole envelope is ejected on a dynam¬ 
ical timescale during CE evolution, we assume that initial 
temperature (To) and density (po) of CE ejecta are equal 
to the temperature and the mass density of stellar matter 
in the giant envelope, respectively. In this work, the stellar 
structure and evol ution are simulated v ia a s tellar evolution 
code developed bv lEggletonl (Il97ll . fl973 ) and lEggleton et al.l 
dl973l) . which has been updated with the latest input p hysics 
over the past three decades (e. g., lHan et all 1 19941 ). The 
giant envelope refers to the region in which the hydrogen 
abundance (by mass) is larger than 0.5. When CE evolution 
occurs, the giant envelope is ejected, shell by shell. Each 
shell expands outwards with no overlapping between shells. 
The To and po of each shell are given by the stellar evolution 
code, and their evolutions are described in the following two 
paragraphs. 

We assume that the CE ejecta is an adiabatically 
expanding perfect gas. The mass density, p, is given by 

p = M ej /(4 7 rR 2 V), (1) 

where M e j, R and V are the mass-ejection rate, the radial 
distance of the ejected matter and the velocity of the ejected 
matter. Some numerical simulations for CE evolution give 
V ~ 100 km s -1 , which is close to the local escape velocity(e. 
g.. lKashi fe Sokerlhoill ). In this work, we assume that 

v * Ksc = y(2) 

where G is the gravitational constant and M is the stellar 
mass. However, considering the conservation of mass during 
CE expansion, we obtain 

M e j = p47rR 2 U = po47rR 2 Vb, (3) 
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where Ro is the initial distance of a given shell of the 
envelope from the stellar center. Here, Vo = sj 2 ^ 1 ■ The 
gas density of the same shell at a distance R from the donor 
center is given by 

p = po(R/Ro)~ 3/2 - (4) 

1 

For an adiabatically expanding perfect gas, pT 1 -Tad = 
constant. Based on the re s ults of model calcu lations in 
iFransson fe Chevalier! Jl989l ') , iKozasa et al.l (Il989l 'l adopted 
an adiabatic index 7 a d as 1.25 for the early stage of SN 
explosion. In this work, we also take 7 a< j as 1.25. Therefore, 
the temperature of the shell at R from the donor center is 
calculated by 

/ p \ -°- 375 

T=To U) • (5) 

Eqs. (1)—(5) determine the evolution of mass density and 
temperature of every shell in CE ejecta. 


2.3 Chemical Abundances of CE Ejecta 


The chemical abundances of CE ejecta are determined by the 
donor’s envelope. For a single star, three dredge-up processes 
and hot bottom burning (the latter only in a star with 
initial mass higher than 4Msun) may change the chemical 
abundances of the stellar envelope. During each dredge-up 
event the base of the convective envelope extends inwards to 
shells previously involved in nuclear burning and nuclearly 
processed ma tter is transported to the externa l zones of 
the star(e. g.. llben fe RenziTilll9831 : lifer widl2005l) . The first 
dredge-up occurs during the first ascent of the giant branch 
(FGB). Sufficiently massive stars (Mi > 4 Mq) undergo the 
second dredg e-up after the end of He-core burning(e. g., 
fHerwig j ' 20051 ). The third dredge-up occurs when the stars 
evolve to thermal pulse AGB (TPAGB) phase. If the mass 
of the hydrogen envelope above the He-exhausted core is 
large enough and the core mass is above ~ 0.8 Mq, the 
hydrogen burning shell can extend into the bottom of the 
convective region, which is called as hot bottom burning. 
Reviews about dredge-ups and hot bottom burn ing were 
given bv llben fc Renziiiil dl98Bh and lHerwiizl (120051 1 . 

As discussed in the previous section, we use Eggleton’s 
code to simulate the evolution of a star from the premain 
sequence to the end of the AGB phase. However, this code 
does not simulate the effects of third dredge-up process and 
hot bottom burning on the che mistry of stellar envelope. In 
this work, we use the model of lLii et al.1 (200811 to simulate 
the chemical evolutions of 4 H, 4 He, 12 C, 13 C, 14 N, 15 N, 16 C> , 
17 0, 20 Ne and 22 Ne in the s tellar envelop e. Lii et al.1 (120081 1 
followed the prescriptions bv llzzard et al.1 (120041 1 for the first 
dredge-up during FGB and the second dredge-up during 
early AGB. For the third dredge-up and the hot bottom 
burning du ring t h e TPAGB phase, we u se the TP AGB 
synth es is in iGro eneweee n fe de J ond (Il9 93h. iKarakas et all 


J'2002il . Ilzzard et al. (l2004l l and Maxigo_&Girardi 


All details can be found in Appendix A of Lii et al 


(2007). 

(2008). 


Mass transfer in a binary system can change the 
chemical abundances of the stellar envelope. In binary 
systems, there are two ways to transfer mass: (i) accretion 
from the stellar wind material of a companion star; (ii) Roche 
lobe overflow. After the gainer in a binary system obtains 


mass ( AM) from the donor, the chemical abundances of 
the stellar surface, X g , are given by 


A' 


new 

g 


X° ld x M| nv + A d x AM 

Mf™T"AM 


(6) 


where M env is the envelope mass and X is the chemical 
abundances of its companion star, g and d indicate the 
gainer and the donor in a binary system, respectively. 

Then, we follow the chemical evolutions of X H, 4 He, 
12 C, 13 C, 14 N, 1S N, 16 0, 1 t O, 20 Ne and 22 Ne in the stellar 
envelope of binary system. The abundance ratio of carbon 
to oxygen, C/O = ( 12 C + 13 C)/( 16 0+ 1 ‘ O), determines the 
dust species produced by CE ejecta. O-rich ejecta(C/0 < 1) 
forms silicate dust, while C-rich ejecta (0/0 _> 1) prod uce 
amorp hous carbon and SiC. Following iFerrarotti fe Gail 
(|200a) . we assume that the abundances of other key elements 
(Fe, Si, Mg and S) for the dust formation do not change. 

Following iLri et al.1 (l2008ll . we use a synthetic stellar 
evolution model to calculate the evolution of chemical 
abundances. Figure [T] shows the time evolution of the 
surface C/O of stars with different masses. The evolutions 
of 0/0 mainly depend on the initial masses of the stars: 

(1) All stars undergo the first dredge-up during the FGB 
phase. Its effects mainly are that carbon abundance on the 
stellar surface is reduced by approximately 30% and oxygen 
abundance does not change. As Figure |T] shows, C /Os of 
all stars decrease after the first dredge-up. The second 
dredge-up, occurring in models with mass above 4 Mq , 
hardly changes the 0/O ratio(e. g. . Ilzzard et al.|[2004| j . 

(2) Low-mass stars (M < 1.5Mq) never become carbon 
stars, owing to the small number of thermal pulses 
experienced. 

(3) Stars with initial masses between about 1.5 Mq and 
4.0 Mq increase their carbon abundance over the oxygen 
abundance after a number of thermal pulses. The surface 
C/O soon exceeds unity. 

(4) Intermediate mass stars with initial masses between 
about 4 Mq and 8 Mg convert the dredged-up carbon 
rapidly into 14 N via the CN-cycle due to hot bottom 
burning. However, the oxygen is not affected by this 
process. Therefore, the carbon abundance in the envelopes 
of these stars is much less than the oxygen abundance. 
These stars never become carbon stars. 


2.4 Dust Formation and Growth in CE Ejecta 

With the decrease of CE ejecta’s temperature due to its 
expansion, dust can form and grow. In this work, with the 
help of AGBDUST code (kindly provided by Prof. Gail, 
private communication), we replace the thermodynamic 
structure of the wind from AGB star in the AGBDUST code 
with Eqs. © and ©, and use the model of formation and 
growth of dust grains in the AGBDUST code to simulate 
the formation and growth of dust grains in CE ejecta. 

The model of formation and growth of dust grains in 
the AGBDUST code assumes that the dust grains grow on 
some kinds of seed nuclei. The seed nuclei are formed prior 
to the point where the different dust components start to 
condense. The radii of these seed nuclei are denoted by ao. 
The AGBDUST code assumes ao = 1 nm for the radii of the 
seed nuclei. If the wind is accelerated by radiation pressure, 
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Figure 1. The evolutions of C/Os in the stellar envelopes with 
different initial masses. The dotted lines represent the evolutions 
of C/Os. The empty and solid triangles give the initial C/Os 
and the C/Os after stars undergo the first dredge-up (FDU), 
respectively. The tops of dotted lines represent the maximum 
values of C/Os at the end of AGB phase. 


Based on the degrees of condensation of key-element (fl) in 
dust species i , the DPR in CE ejecta is calculated by 


dM d _ , > y ^ fi 

dt - M ei X k ^ fk, 


( 11 ) 


where X k is the key-element abundance in CE ejecta by 
mass, Ai and Ak are the molecular and atomic weights of 
the dust species i and the key-element fc, respectively. 

Tabic [T] gives a list of the dust species considered in the 
present paper, their formation reactions, the corresponding 
key-elements and the sticking coefficient {ai). The ther¬ 
modynamic quantities (mainly the free energies used to 
calculate t he equilibrium pre s sures) of the reactions are 
taken from ISharp fe Huebneil (1 19901 1. The composition of 
olivine and pyroxene depends on the relative fraction of 
magnesium and iron, indicated, with x and 1 — a: (see Table 
Q}, respectively. The details for calculating the magnesium 
percentages in the olivine an d pyro xene dust grains can 
be found in iGail fe Sedlmavij (| 19991 ). To our knowledge, 
there is not experimentally determined da ta for the sticking 
coeffic ient. We adopt the values of ai in Ferrarotti fe Gaill 
(|2003|) which are mainly based on lHashimotol 1 1990l ) and 
iNagahara fc Ozawal Jl99fil ). 


the final results for the grain radii are nearly independent 
of the particular assumption on ao . Then, the growth of the 
grain size of species i is calculated via the equation 

^ = v 0>i {jr - Jtn, ( 7 ) 

where ai is the radius of dust species i, and Vo,; is the 
volume of the nominal molecule in the solid. Here, Jf 1 is the 
deposition rate per unit time and surface area for a particle 
in rest with respect to the gas phase. It is given by 

Jf = aimvth.i, (8) 


where a; denotes the sticking coefficient, m gives the 
particle density in the gas phase and u t h,; = is 

the mean thermal velocity of the particles. Here, k is the 
Boltzmann constant and T is the temperature. The jf ec is 
decomposition rate in a thermodynamic equilibrium state, 
which is given by 

Pi 


J dec 
i — 


kT’ 


(9) 


where p; is the partial pressure of the molecules of dust 
species i in chemical equilibrium between the gas phase and 
the solid phase. 

The amount of dust formed is indicated by the frac¬ 
tion of gas co ndensed into solid grai ns, i.e. the degree 
of condensation llGail fc Sedlmavii 119991) . For each species, 
the AGBDUST code selects a key-element k. Its degree of 
condensation is expressed by, fl , the fraction of the key- 
element in dust species i, which is given by 


i _ 47r[(a;) 3 - (a 0 ) 3 ] n d> ; 
lk ~ 3Vq,; t k Nn ’ 


( 10 ) 


where n d ,; is the number density of seed nuclei of the 
different dust species, e k is the key-element abundance by 
number relative to H abundance, and TVh is the number 
density of hydrogen nuclei. In the AGBDUST code, w d .; = 
3 x 10~ 13 JVh for all dust species dFerrarotti fe Gail 120031) . 


3 DUST YIELDS FROM CE EJECTA 

Based on the model of dust formation and growth in CE 
ejecta described in the last section, dust yields not only 
depend on the parameters listed in Table [l] but also on 
the initial temperature (To) and density (po), the chemical 
abundances (mainly C/O) and the mass of CE ejecta. 

In this section, we take Z = 0.008 as an example 
to calculate the dust yields from CE ejecta. The initial 
abundances of all heavy elements are approximately equal 
to 0.4 times those of the Sun. The ele ment abu nd ances on 
the su rface of the Sun are taken from lAnders fc Grevessd 
d!989P . 

3.1 Typical Examples 

As discussed in (J21 CE evolution involves a giant or giant-like 
envelope. This giant star may evolve in the Hertzsprung gap 
(HZ), FGB or AGB. The envelopes at different evolutionary 
stages have different structures (mass density, temperature 
and chemical abundances), which may affect the dust 
formation and growth in the CE ejecta, fn order to discuss 
these effects, we select two groups of typical examples: one 
is that CE ejecta originate from an envelope of a star with 
initial mass of 1 Mq at the beginning of HZ, the beginning of 
FGB, the end of FGB, the beginning of AGB and TPAGB, 
respectively; the other is that CE ejecta is produced by an 
envelope of a star with initial mass of 3 Mq at the begin 
of HZ, the begin of FGB, the end of FGB, the begin of 
AGB, and C/O = 2 at TPAGB due to the third dredge-up, 
respectively. 

Figure 0 shows the initial temperature and density 
profiles of the giant envelopes as a function of relative mass 
coordinate for the two groups of typical examples. The 
external regions of the stars expand during the evolution, 
with the consequent decrease of temperature and density. 
The mass of the envelope diminishes, owing to mass loss. 
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Table 1 . Dust species considered in the present paper, their formation reactions, the corresponding key-elements, and the sticking 
coefficient (ft,). 


Dust Species 


Formation Reactions 

Key-Elements 

OLi 

Olivine 

2zMg + 2(1 

— x)Fe + SiO + 3 H 2 O —y Mg 2 a;Fe 2 (i_ a ;)Si 04 -f- 3 H 2 

Si 

0.15 

Pyroxene 

xMg+ (1 

— :r)Fe SiO + 2 H 2 O —y Mga;Fe^—x)Si 03 + 2 H 2 

Si 

0.15 

Quartz 


SiO + H 2 0 —>• Si0 2 (s) + H 2 

Si 

0.1 

Silicon Carbide 


2Si + C2H2 ->• 2SiC + H 2 

Si 

0.3 

Carbon 


C ->■ C(s) 

C 

0.3 

Iron 


Fe —► Fe(s) 

Fe 

0.8 


As shown in Figure[5] at the beginning of CE expansion 
the temperatures in the envelope are too large (higher than 
~ 10 4 K) to allow dust formation. During the expansion the 
CE ejecta reach regions 10 13 —10 16 cm away from the donor, 
where the temperatures are below 2000K, thus rendering 
possi ble condensation of gas molecules into dust dLii et ali 
l2013l l. By means of the AGBDUST code we simulate dust 
formation and growth within the expanding envelope. The 
computation is run for each shell, of initial temperature To, 
into which the envelope is divided. Figures [3] and 0] show the 
condensation degrees of key elements and the mass of dust 
formed in the expanding shells of the envelope of two stars 
of mass 1 and 3 solar masses. 

The dust formation and growth in CE ejecta has its 
own features: 

(i) High condensation degrees. The main reason is that the 
dust-forming zones in the CE ejecta (T <~ 1300 K) have 
very high mass density (~ 3.3 x 10” n g cm” 3 ). Based on 
Eqs.([7]) and l[8]l. the higher is the mass density, the higher 
is the growth rate of grain size, which results in a high 
condensate degree. For AGB wind with a high mass-loss 
rate of 1 x 10~ 5 Mpiyr~ 1 and a wind structure based on 
iFerrarotti fe Gail (l2006l l . AGBDUST code shows that the 
temperature and the mass density of dust-forming zones are 
~ 1.3 x 10 3 K and 5.0 x 10” 14 g cm” 3 , respectively. About 
32% of Si atoms in AGB wind condensate into silicates, 
and about 4% of Fe atoms condensate into iron grains. 
However, in our simulations, all of Si and Fe atoms in CE 
ejeta condensate into silicate (/ol + /qu + /py ~ 1) and 
iron grains (fm ~ 1), respectively. 

(ii) Silicate grains are mainly made up of pure forsterite 
(Mg 2 SiC> 4 ) and quartz grains. At the region of dust 
formation and growth, compared with AGB stellar wind, 
CE ejecta has higher temperature and mass density. Based 
on the sta bility limits of some gra ins and molecules (See 
Fig. 2 in I Gail fe Sedlmavil (j 19991 4 4. forsterite, iron and 
quartz grains firstly form. Most of Fe atoms condensate into 
iron grains. Therefore, olivine-type (Mg 2 a: Fe 2 (i_ a; )Si 04 ) 
grains are mainly made up of pure forsterite (Mg 2 Si 04 ) 
grains. 


As Figures [3] and [4] show, almost all of Fe atoms in CE 
ejecta condensate into iron grains. In a CE ejecta ( C/O < 1), 
most of Si element condensates into silicate grains; In CE 
ejecta ( C/O > 1), a part of C atoms and most of O 
atoms combine CO molecule, while most of left C atoms 
condensate into carbon grains. Therefore, CE ejecta can 
efficiently produce dust grains, and the dust mass mainly 
depends on the mass of CE ejecta. 


3.2 Comparison with AGB Stars 


AGB stars, as important stellar sources of dust, have been 
investigated by many independen t groups. Havin g used 
synthetic stellar evolution models, IZhukovska et all d2008l l 
computed dust condensation and estimated the dust yields 
for M, S, and C-type AGB stars of different ages and 
metal licities. Based on these models, IZhukovska fe Henning! 
(|2013ih esti mated dust input from AGB stars in the LMC. 
However, as I Schneider et al] d2014l 4 pointed out, their major 
drawback was that they used simple evolutionary code 
based on analytic approximations to describe complex phys¬ 
ical processes, primarily third dredge-up and hot bottom 
burning. These processes directly affect the chemical abun¬ 
dances and physical conditions in AGB envelope, thus 
they determine the dust yields. Ventura and collaborators 
used the ATON code to produce accurate modelling of 
the AGB phase and computed the mass and com p ositio n 
of dust produced b y AGB stars ( VentumetalJ d2012li ; 
IVentura et al.l (20121): Di Criscienzo et al.l 1 201314 4. Follow¬ 


ing 


Schneider et al. (2014), we call this code as old ATON. 


IZhukovska fe Hennine ( 20131 4 found that, compared with 
observations, the old ATON yields underestimated the dust 
production ra t e from C-rich AGB stars. For this reason, 
IVentura et alJ (120141 4 calculated new AGB evolutionary se¬ 
quenc es with a deeper third dredge-up, and lSchneider et al.l 
(|2014l 4 provided new dust yields, to which we refer to as 
new ATON. However, there are some differen c es of dust 
mass and composition among IZhukovska et al.l (l2008l 4. old 
and new ATON. Det ailed reasons had been d iscussed by 
ISchneider et al.l d2014l4 and IVentura et all (120141 4 . 

As discussed in next section, based on the observational 
data, several investigations estimated the DPRs from AGB 
stars in the LMC, and their results have large scatter. 
By comparing the observed and pr edicte d DPRs of the 
Magellanic Clouds, ISchneider et al] (|2014| j found that old 
ATON models are better suited for Z <= 0.004, while use 
of the new ATON models is recommended for Z = 0.008. In 
this work, in order to compare the dust yields produced by 
CE ejecta with those by AGB stars, we use the above three 
models. 

In Figure [5J we compare the dust masses of different 
compositions from CE ejecta formed at di fferent evolution¬ 
ary st ages with those from AGB stars bv IZhukovska et al.l 
(boosh . old and new ATON. Due to the high condensation 
degrees in almost the whole envelope, CE ejecta produces 
more efficiently silicate and iron grains. In our model, the 
stars with initial masses between ~ 1.5 and 4.0 Mq evolve 
to carbon stars and produce carbon and SiC grains, similar 
to the results of old and new ATON. On the other hand, 
having used a model for TPAGB with efficient third dredge- 
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Figure 2. The initial temperature and density profiles of the envelopes of giants (at different evolutionary stages as indicated) as 
a function of relative mass coordinate. ’HZb’, ’FGBb’, ’FGBe’, ’EAGB’ and ’TPAGBb’ mean that star just evolves at the begin of 
Hertzsprung gap, the begin of the FGB, the end of the FGB, the begin of the AGB and the begin of the thermally pulsing AGB 
(TPAGB), respectively. ’TPAGB2’ means that C/O of stellar envelope increases to 2.0 because of the third dredge-up. The left and right 
panels are for the star with mass of 1 and 3 Mq at zero-age main sequence, respectively. 




M r (M & ) 


Figure 3. Condensation degrees of key-elements and dust yields of different dust species as a function of relative mass coordinate at 
different evolutionary stages for the star with mass of 1 Mq when stellar envelope is ejected into a region where dust can form and 
grow(Details in text). The ’/ol\ ’/qu’ and 7py’ present the condensation degrees of Si element in olivine-type, pyroxene-type and 
quartz grains, respectively. The ! /lIt' means the condensation degree of Fe element in iron grains. Correspondingly, the ’A/<j’s present 
the dust yields of different dust species. 


up and w eak hot bottom b u rning , the stars between ~ 1 and 
8 Mq in IZhukovska et all ( 2008 ) can become carbon stars 
with high C/O at the end of AGB stars. Therefore, their 
models produced more carbon grains and less silicates. 


3.3 Comparison with Supernova 

SNe II are considered as another important stellar sources 
of dust. Due to their short lifetimes and large production 
of metals, they are regarded as the most likely s ources 
of dust in high red-sh ift galaxies and quasars (e.g., iDwekl 
1 19981 : lOall et all 1201 ll ). Unfortunately, no matter observa¬ 
tions or theoretical models, there is no reliable result to 
definitively calculate which species of dust formed and in 
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Figure 4. Similar to Figure 0 but for the star with mass of 3 Mq. The ’/ca’ and ’/sic’ present the condensation degrees of carbon 
and Si elements in carbon and silicon carbide grains, respectively. 
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Figure 5. Comparison of the dust yields fr om CE ejecta formed at different evo lutionary stages with AGB stellar winds . Z2008, old and 
new ATON mean that the data come from lZhukovska et alj d2008h . old ATON dVentura et al.l (1 201 2| ) ) and new ATON llSchneider et all 
^2014j V). respectively. 


what quantities by SNe II. Several theoretical models for 
dust formation in the ejecta of SNe II have been developed 


fe.g.. ITodini fe FerraralfeoQll: ; NozaMmet_aL 20031'). Based on 
theory, [ 


standard nucleation theory. iTodini fc FerrarJ ~( 200 il l calcu¬ 


lated dust yields in the ejecta of SNe II whose progenitors 
have different metallicities (Z = 0, Z = 10 -4 , Z = 1CP 2 and 
Z = 0.02). By a linear interpolation between Z = 10~ 4 


and Z = 10~ , we can obtain dust yields from SNe II 
whose progenitors have Z = 0.008. Ev ery S N can produce a 
mass of dust of up to 1 Mq. Howe ver, Bianchi_J^_Schneidei| 
d2007t) reconsidered the model of ITodini fc Ferraral ( 20oiF ) 
and included dust destruction. They found that only about 
2 — 20 percent of the initial dust mass can survive the 
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Figure 6. Comparison of the dust yields from CE ejecta 
formed by massive stars at the begin of red supergiant (RSG). 
The dust yields f rom SNe II without reverse shock come from 
iTodini &; Ferrari (|2Q0ll j. The dust yields from SNe II with reverse 
shock mean that only about 2 percent of dust mass survive after 
reverse shock. 

passage of the reverse shock, depending on the density of 
the surrounding interstellar medium. 

In a binary system, if a massive star fills its Roche 
lobe when it evolves into the red supergiant (RSG) state, 
its envelope can also be ejected as CE ejecta. We compare 
in Figure[6]the dust yields in CE ejecta with those from SNe 
II. In the ejecta of SNe II, there is a large quantity of metal 
elements. In the model of ITodini fc Fe rraral (1200lh . most of 
them condensate into dust grains. Therefore, compared with 
the ejecta of SNe II without dust destruction, dust produced 
b y CE ejecta from massiv e star is negligible. However, as 
iBianchi fc Schneider! d2007l l showed, most of dust produced 
by the ejecta of SNe II are soon destroyed by reverse shock. 
If only about 2 percent of the initial dust mass produced by 
the ejecta of SNe II can survive after the reverse shock, dust 
produced by CE ejecta from massive star is dominant. 


4 DUST BUDGET IN THE LMC 

As shown in the last section, AGB stars, the ejecta from 
SNe II and CE are potential sources of dust in the LMC. 
In order to calculate their contributions to total dust in the 
LMC, we use a method of population synthesis. 

4.1 Basic Parameters of Population Synthesis 

For a method of population synthesis, the main input pa¬ 
rameters are:(i)the initial mass function (IMF); (ii)the mass- 
ratio distribution for binaries;(iii)the eccentricity distribu¬ 
tion for binaries; (iv)the distribution of orbital separations 
for binaries. We use a Salpeter initial mass-function for 
single stars and primary components in binaries: 

4>{m) = Cm~ 2 ' 35 , (12) 

with normalisation number C = 0.06 within the mass 
integral range between 0.1 and 100 Mq. For binary systems, 
we must know the distribution of the ratios of components. 


Based on iKraicheva et al.l (Il989l l and iGoldberg fc Mazebl 
d 1994ll . a flat distribution of initial mass ratios of compo¬ 
nents is assumed in this work. Meanwhile we assume that 
all binaries have initially circular orbits. In our work, one 
of the most important input parameters is the distribution 
of initial binary separations which directly determine how 
many binaries undergo Roche-lobe overflow. However, to our 
knowledge, the re is no a n y conf i rmed observational evidence 
to describe it. lEldridee et alj (I2008T ) assumed that a flat 
distri bution over 1 < log ao/ Rq < 4 (also see iTout et alJ 
l2014ll . where ao is the initial binary separation. In this work, 
we adopt the above distribution. 

4.2 Star Formation Rate 

The dust quantities in the LMC depend on the star forma¬ 
tion rate (SFR). According to the Magellanic Clouds Pho¬ 
tometric Survey including 20 mil lion LMC stars an d using 
the StarFISH analysis software, [Harris fc Zaritskvl (1200911 
reconstructed the SFR of the LMC during different ages. 
[Harris fc Zaritskvl (120091 ') found that the optimal solution in¬ 
cludes the four metallicities Z = 0.001,0.0025,0.004,0.008. 
The separate contributions of the four metallicities during 
different ages are shown in Figure [7] Since we are interested 
in the dust quantities of the whole LMC, we average over 
SFRs of individual regions for each age bin. 

4.3 Initial Binary Fraction in the LMC 

In our work, dust formation involves binary systems. There¬ 
fore, the binary fraction in all stars of the LMC is very 
important. There is no obs ervational evidence to he lp us in 
adopting a binary fraction. lHarris fc Zaritskvl (|2009l ) simply 
assumed that the binary fraction in the LMC is 50%. 
However, it does not give solid information for the initial 
binary fraction. According to the theoretical model of binary 
evolution, binary systems can become single stars if they 
undergo binary merger or they are disrupted by supernovae. 
For example, for binary systems whose initial conditions are 
given in 34.11 about 38% of binaries do not undergo Roche 
lobe mass transfer, and their evolutions are similar to single 
stars. About 62% of binaries undergo Roche lobe overflow 
at lease once, in which about 73% undergo CE evolution 
and can eject the whole CE at lease once, while about 27% 
can not eject the CE and usually merge into single stars. In 
short, about 54% of these binary systems may be observed 
as single stars. 

In this work, in order to check the effects of initial binary 
fraction on the dust budget in the LMC, we use two different 
initial binary fractions: one is an initial binary fraction of 
50% (a half of SFR in the LMC is for binary systems and 
another half is for single stars), the other is that the initial 
binary fraction is 100%, that is, all stars are born in binary 
systems. 

4.4 Calculation of Dust Production Rates 

Using IMF, SFR and dust yields discussed in the previous 
sections, we can theoretically calculate DPRs from the above 
three sources in the LMC. _ 

For a single star, following IZhukovska fc Hennind 
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Figure 7. Star fo r matio n history in the LMC with the age from 
I Harris &: Zaritskvl d2009T) . The best-fit SFR as a function of age is 
showed with a solid line, and the uncertainty on the fit is showed 
as an oblique-line area. Dashed, dash-doted, doted, dash-dot-dot- 
doted lines present the contributions of Z = 0.001, 0.0025, 0.004 
and 0.008, respectively. 


(| 20131 ) and [Schneider et al.l ([20141 ), the DPRs of the type j 
dust grains from AGB stars at time t can be estimated by 




/*n 

J m 


<P(m] 


SFR(t — r(m, Z)) 


m d (m, Z)dm, 


(13) 


where 4>(m) is IMF, SFR(t — r(m,Z)) is SFR, T(m,Z) 
is the lifetime of a star with mass m and metallicity Z 
when it leaves AGB stage, and m d (m,Z) is the mass of 
dust. The lower integration limit, itil, is determined by a 
stellar mass with the lifetime which satisfies r(mi) = t, and 
the upper mass-limit mu = 8M 0 is determined by stars 
becoming AGB stars. The m d (m,Z) can be obtained by 
the linear interpolation i n a given grid of dust yields from 
IZhukovska et al] <2008l ) or lSchneider et ahl <20141 ). Similarly, 
the DPRs of the type j dust grains from SNe II at time 
t can also be calculated by Eq. Gi, but the upper mass- 
limit rri\j = 40M 0 and the m d (m, Z) can be obtained from 
iTodini &; Ferraral (120011 1. 

For a binary system in which the initial masses of 
primary and secondary are mi and m 2 and the initial binary 
separation is ao, the DPR of the type j dust grains from CE 
ejecta at time t can be estimated by 


Mi = 


iz L v ^( m i) sz1 ^( m2 ) /io° 

SFR(t-T(m 1 ,m 2 ,a 0 )) ( mi; m2; ao )dm 1 dm2da 0 


where r(mi, m 2 , a 0 ) is the time of a binary system occurring 
CE evolution, and mJ(mi,m 2 ,ao) is dust mass produced 
by CE ejecta. In order to calculate m d (mi, m 2 , ao), for 
the four metallicities (Z = 0.001,0.0025,0.004,0.008), we 
calculate dust yields which are for different stellar masses 
(m = 1.0,1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0, 5.5, 6.0, 8.OM 0 ) in 
different evolutional phases(HZ, the beginning of FGB, the 
end of FGB, the beginning of AGB and TPAGB). We do not 
consider massive stars due to mass function. When a binary 
system undergoes CE evolution, by the linear interpolation 
of the donor’ mass and age in the dust yields, we obtain 
m d(mi, m2, 00). 


4.5 Dust Production Rates 


Based on observational data of a number of surveys, DPRs 
from C-ric h and O-rich AGB stars in the LMC had been es ¬ 
timated by Matsuura et a l.l <2009l) .ISrinivasan et alJ (2009), 


iBover et al 1 20121 ) and Riebel et alJ 12012 tl . Unfortunately, 
as Figure [5] shows, there is a large scatter between the 
above observationally est imated DPRs. DPRs f rom C-rich 
AGB stars estimated by iMatsuura et al.l <2009l ) are larger 
than 4.3 x 10 _5 M 0 yr -1 (up to 1. 0 x 1O~ 4 M 0 yr -1 ), while 
they are 2.4 x 10~ 6 M 0 yr -1 in ISrini vasan et ahl 1 20091) . 


8.7x10“ 


% y f 


yr 1 in iRiebel et al 


Bover et al. <20121) . 5.2—5.7x 10" 6 M m 
( 201 2l ). respectively. Similar scatter 
appears in DPRs from O-rich AGB stars. Differences of the 
DPRs observationally estimated mainly result from various 
approaches and source classif ications (Details can be seen in 
IZhukovska fe Hennind (|2013l )). 

Similarly, there is also a large difference between the 
theoretica lly estimated DP R s. Usin g the dust yields of AGB 
stars in [Zhukovs ka et al] <20081 ). IZhukovska &; Hennind 
(|2013i) predicted that the DPRs from C-rich and O- 
rich AGB stars in the LMC are 5.7xlO~ 5 M 0 yr -1 and 
1.3 x10~ 6 Mq yr -1 , respectively. Ho wever, using the dust 
yields of AGB stars from old ATON. IZhukovska fe Hennind 
<2013l) calculated that they were 6.3xlO~ 6 A/ 0 yr -1 and 
1.1 x10~ 5 Mq yr -1 , respectively. The uncertainty introduced 
by different models for AGB stars is up to ab out 1 0. Usin g 
the dust yields of A GB stars from Zhukov ska et all (2008), 
old and new ATON. I Schneider et al.1 ( 20141 ) also predicted 
the DPRs in the LMC. They found that the DPR from 
C-rich AG B stars calculated by the dust yields of AGB stars 
from IZhukovska et al] (|2008l ) is close to that calculated by 
new ATON, while the DPR from O-rich AGB stars by new 
ATON is close to that by old ATON. The main reason 
behind these differences is in the treatment of convective 
borders and in the efficiency of the convective models 
used, that affect the extent of the third dredge-up and 
th e strength of ho t bottom burning (Details can be seen 
in ISchneider et al.l <2014l )). However, to our knowledge, 
all theoretical estimates did not consider the effect of 
binary interaction on the DPRs in the LMC. In this work 
we compar e the D PRs of A GB s tars in the LMC with 
those from IZhukovska et al] <2008l) . old and new ATON. 
Simultaneously, we also estimates the DPRs of CE ejecta 
in the LMC. As Figure [8] shows, the initial binary fraction 
changing from 50% to 100% introduces an uncertainty 
of a factor of about 1.5 — 2 to the DPRs from AGB 
stars and from CE ejecta, respectively. These uncertainties 
are much smaller than these differences between different 
observational or theoretical estimates. 


Due to large scatters between observationally estimated 
or theoretically predicted DPRs from AGB stars, we can not 
determine whether the DPRs from CE ejecta are significant. 
However, we find that CE ejecta can efficiently produce 
silicate and iron grains, while carbon grains are negligible. 
The main reason is that in our models most of CE events 
occur before the third dredge-up works efficiently. 

I 11 order to compare the contributions of different 
sources to the DPRs in the LMC, Figure [9] shows the 
total DPRs from CE ejecta, AGB star and SNe II with 
the age of the LMC by using the best-fit star formation 
histories. If we do not consider the destruction by reverse 
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shock, the majority of dust grains in the LMC originate 
from SNe II. However, if only 2% of dust grains can survive 
after reverse shock, com pared with AGB stars calculated by 
IZhukovska et al] (l2008ll and new ATON, the contribution 
of SNe II to the total DPRs is negligible. In the simulation 
with an initial binary fraction of 50%, the DPR from CE 
ejecta is close to that calculat ed by old ATON. but it is 
only 1/4 of those calculated bv IZhukovska et al.l (1200811 and 
new ATON. However, it is about 4 times of the former and 
about 3/4 times of the latter in the simulation with an 
initial binary fraction of 100%. One should take note that 
the assumption of 100% binary fraction is realistic. Here, 
we only want to estimate an upper limit of the DPR from 
CE ejecta. More imp o rtantl y, most of dust grains calculated 
by IZhukovska et al] (l2008h and new ATON are carbon. 
Formation of C-rich AGB stars deeply depend on a complex 
phy sical process: the t hird d redge-up. In new ATON model 
and IZhukovska et al] (l2008l l. the third dredge-up is more 
efficient, and C-rich AGB stars form more easily. However, 
in old ATON, it is difficult for a star to evolve into C-rich 
AGB star, and the DPR from C-rich AGB stars in the LMC 
is negligible dZhukovska et al1l2008l : ISchneider et al-lliohlf . 

In short, although the DPRs in the LMC estimated by 
observational data or theoretical models are very uncertain, 
we consider that the contributions of CE ejecta must be 
included when the total DPRs in the LMC are calculated. 


4.6 Total Dust Budget 

The total dust mass, Mdust, is very important for under¬ 
standing the dust formation, evolution and destruction in 
the LMC. Therefore, there are many literatures to estimate 
M dust in the LMC by using some observational datafe . 
g., iBot et ai]l20ld : ISkibba et al.ll2012l : I Gordon et al. 1 12014) . 
However, the total dust masses estimated are different: 
Using sub- millimeter ex c ess of t he LMC and fittin g the dust 
models of iDraine fc H (l2007h . IBot et all (j20ld l obtained 
Mdust of 3.6 x 10 6 Mq; A ccording to Herschel HERITAGE 
data. ISkibba et all d2012il predicted Mdust ~ 1.1 X 10 6 Mq; 
Base d on the HE RIT AGE Herschel Key Project photometric 
data. iGordon et ahl (|2014l j obtained Mdust of 7.3 x 10 5 ALq, 
which is a factor of about 5 lower than previous results. 
These differences maybe come from different dust models, 
the fitting techni ques, the wavelength range, and the factor 
of environments d Gordon et al.|[2014l h As Figure flOl shows, 
the total dust mass produced by AGB stars or CE ejecta 
is much lower than the above observational estimates. Of 
course, we do not know whether dust produced by SNe II 
is dominant because theoretical models and observational 
estimates can hardly give any certain value for the dust 
yields of SNe II. 

Furthermore, dust grains can be destroyed by the blast 
wave produced by SNe. Therefore, their existence should 
have time scale. According to the properties of the dust 
material and the occurrenc e rate of SNe in the LMC, 
IZhukovska fe Henningj d2013T ) estimated the time scale of 
dust destruction. They found that the time scale for silicate 
grains is about 0.8 Gyr, and it is about 1.1 Gyr for 
carbonaceous grains. This means that about 10% of dust 
masses showed in Figure [TO] can survive if we consider the 
dust destruction. It is possible that AGB stars and CE ejecta 


are not m ain dust sour ces and dust can grow in the ISM of 
the LMC dDrainell2009l h 

As discussed in I ntrod uction , according to the estimates 
of IWeingartner fc Drainel (l200lT ) , we indirectly guess that 
the mass ratio of carbon to silicate grains in the LMC is 
about 1:4. However, to our knowledge, there is no direct 
observations to give the compositions of different dust 
species in the LMC. It is very important for us to understand 
the dust origin. We theoretically predict the time evolution 
of accumulated masses of different dust species from CE 
eje cta and AGB stars in F i gure Til Obviously, in the model 
of IZhukovska fe Henning! ( 2013 ), most of dust grains are 
carbon and others are negligible. In the model of new 
ATON, the mass fractions of carbon and silicate grains 
are about 82% and 13%, respectively, and the left are 
mainly SiC grains. In the model of old ATON, the mass 
fractions of carbon and silicate grains are about 55% and 
35%, respectively. The remaining are SiC and iron grains. 
However, in CE ejecta, about 50% of dust grains are silicate, 
20 % of them are iron and the mass fractions of carbon 
and SiC grains are about 17% and 13%, respectively. Not 
like AGB stars which mainly form C-rich grains, CE ejecta 
chiefly produces O-rich dust grains. 

Therefore, with the further observations, we believe that 
the compositions of different dust species in the LMC will be 
found soon. Then, we can estimate which between AGB stars 
and CE ejecta is more important dust source. Of course, we 
also should consider the effect of SNe II and dust growth in 
the ISM. 


5 CONCLUSIONS 

In this work, we calculate the contributions of dust produced 
by AGB stars, CE ejecta and SNe II to the total dust budget 
in the LMC. Because different observations still have large 
scatter and theoretical models greatly depend on complex 
physical processes, we can not determine which is the main 
dust source in the LMC. However, more than 50% of stars 
in the LMC are in binary systems. Many of them may have 
or will undergo CE evolutions. According to our model, 
CE ejecta can efficiently produce dust. Therefore, we must 
consider the contributions of dust produced by CE ejecta on 
estimating the total dust budget. 

In our simulations, the DPRs of AGB stars in the LMC 
are between about 2.5x 10 _5 MQyr _1 and 4.Ox lO~ 6 M 0 yr _1 . 
The uncertainty mainly results from different models for 
the dust yields of AGB stars. These results are within 
the large scatter of several observational estimates. The 
DPRs of CE ejecta are about 6.3 x lO^MgyU 1 (The 
initial binary fraction is 50%). Compared the theoretically 
predicted DPRs of AGB stars, the DPRs of CE ejecta are 
significant or dominated. The contribution of SNe II is very 
uncertain. Compared with SNe II without reverse shock, the 
DPRs of AGB stars and CE ejecta are negligible. However, 
if only 2 % of dust grains produced by SNe II can survive 
after reverse shock, the contribution of SNe II is very small. 

The total dust mass of the LMC estimated by obser¬ 
vations is between 7.3 x 10 s M@ and 3.6 x 10 6 Mq. In our 
work, the total dust mass produced by AGB stars in different 
models are between 2.8 x 10 4 Mq and 3.2 x 10 5 Mg, and 
those produced by CE ejecta are about 6.3 X 10 4 Mq. In 
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Figure 8. The dust production rates (DPRs) of different dust species from CE ejec ta and AGB stars in the LMC as a function of time 
for metallicity-dependent star formation histories given bv lHarris & Zaritskvl (|20091 ). Left and right panels show that the initial binary 
fractions are 50% and 100%, respectively. The solid, dashed, dash-d oted and dot-dot-dot-dashed lines show the predicted DPRs from CE 
ejecta, AGB stars using the dust yields of lZhukovska et al l (|2008l ). new and old ATON, respectively. The shaded regions around lines 
re present the uncertainty on the best-fit star formation histories. The DP Rs of O-rich and C- rich AGB stars observationally estimated 
bv lMatsuura et al.1 <[20091 ). ferinivasan et al . ( 20091). bover et al.1 (|2012h and lRiebel et all (|2012h are represented by triples, stars, squares, 
and circles, respectively. 




Figure 9. The total DPRs from CE ejecta, AGB stars and SNe II with the age of the LMC by using the best-fit star formation histories. 
Left and right panels show that the initial binary fractions are 50% and 100%, respectively. Th e solid, dashed, dash-do ted and dot- 
dot-dot-d ashed lines show the p redicted DPRs from CE ejecta, AGB stars using the dust yields of lZhukovska et al.l l2008|) and the new 
ATON in lSchneider et abl ll2014h . respectively. Thin doted l ines represent the DPRs from SNe II using the dust yields of lTodini &; Ferraral 
[1 2001 1 ) . while thick doted lines use 2% of the dust yields of lTodini fe Ferraral [1 2001 1 ) . 


our calculations, the total dust masses produced by SNe II 
are vary uncertain, and they change from 2.8 x 10 3 Mq to 
2.8 x 10 6 Mq. If the destruction of blast wave from SNe is 
considered, compared with the observational estimates, the 
total dust masses produced by AGB stars, CE ejecta and 
SNe II are negligible. Therefore, on estimating the total dust 
mass, the growth of dust in the ISM must be considered. 

Not only observational estimates but also theoretical 
predictions found that AGB stars in the LMC mainly 
produce carbon grains. Based on our simulations, most of 
dust grains forming in CE ejecta are silicate and iron grains. 
Therefore, with the further observations, the compositions 


of different dust species will be found, which may be 
very important for understanding the dust origins and the 
contributions of dust produced by AGB stars and CE ejecta. 
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Figure 10. Similar with Figure |9| but for time evolution o f total dust masses from CE ejecta, AG B s tars and SNe I I . The observational 
estimates presented by circle, square and star are taken from lBolT et alJ ((201Qh . lSkibba et al.1 (|2012l) and I Gordon et ahl (|2014h . respectively. 




Figure 11. Similar with Figure[8] but for time evolution of accumulated masses of different dust species from CE ejecta and AGB stars. 
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